Anticorrelation between Ion Acceleration and Nonlinear Coherent 
Structures from Laser-Underdense Plasma Interaction 
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In laser-plasma experiments, we observed that ion acceleration from the Coulomb explosion of 
the plasma channel bored by the laser, is prevented when multiple plasma instabilities such as 
fllamentation and hosing, and nonlinear coherent structures (vortices/post-solitons) appear in the 
wake of an ultrashort laser pulse. The tailoring of the longitudinal plasma density ramp allows 
us to control the onset of these insabilities. We deduced that the laser pulse is depleted into 
these structures in our conditions, when a plasma at about 10% of the critical density exhibits 
a gradient on the order of 250 jim. (gaussian fit), thus hindering the acceleration. A promising 
experimental setup with a long pulse is demonstrated enabling the excitation of an isolated coherent 
structure for polarimetric measurements and, in further perspectives, parametric studies of ion 
plasma acceleration efficiency 
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Laser-underdense plasma acceleration rests on the in- 
teraction of an intense relativistic pulse (IX 2 > 10 18 
Wcm _2 /im 2 ) with a plasma in a large range of electronic 
density (n e G [10 18 ,10 20 ] cm" 3 ) [lffi]. In particular, the 
Coulomb explosion of a positively charged channel, cre- 
ated by the laser as it propagates in the plasma, is an 
efficient way to radially accelerate ions [U[5]. However, 
this acceleration can be accompanied by the copious gen- 
eration of non-linear plasma structures that significantly 
deplete the pulse energy, theoretically up to 30 — 40% of 
the laser energy [6UTU] . These structures aroused the in- 
terest of numerous theoretical and experimental studies 
[6j [HJ [TTHT5] . An underdense plasma with n e ~ tens of 
percent of n c , n c being the critical density, is particularly 
suitable for the excitation of coherent structures like elec- 
tromagnetic solitons and electron vortices, as predicted 
in Refs. [6j [8] [TO] [11] and observed with proton radio- 
graphy diagnostics in Refs. p~2j [13j [16]. In Ref. [6] a 
3D electromagnetic cavity structure in the wake of the 
pulse, termed transverse magnetic-soliton, is identified 
and slowly expands into a post-soliton, as ions within the 
cavity explode to energies of tens of keV. This example 
illustrated a novel ion acceleration scheme from under- 
dense plasma, different from the channel Coulomb explo- 
sion PQ, Target Normal Sheath Acceleration-like process 
[3], or longitudinal magnetic vortex acceleration [T7] . 

In this letter, we present for the first time an "inverse 
correlation" (mutual exclusion) between the ion acceler- 
ation by radial Coulomb explosion and the generation 
of nonlinear coherent structures in the wake of the laser 
pulse. Our results tend to prove that the efficiency of 
the laser ion acceleration directly from a gas jet is highly 
dependent upon the density gradient profile. In our case, 
though the density ramp should contribute to the in- 
crease of the laser intensity by self-focusing of the pulse 
P~8j [19], it induces fllamentation [19 , hosing instabili- 



ties [20] and ultimately generation of nonlinear coherent 
structures, the overall turning out to be detrimental for 
the acceleration. 

The experiment was conducted at the Laboratoire 
d'Optique Appliquee with the Ti: Sapphire laser "Salle 
Jaune", which delivers an ultrashort (duration r = 30 
fs), linearly-polarized (along z, see FIG. [TJ pulse with 
El ~810 mJ on target [21]. The pulse was focused with 
an //10 off- axis parabolic mirror to a 25-/im (1/e 2 ) fo- 
cal spot onto a 700 jam (FWHM) supersonic He jet at 
~300 /im above the nozzle exit. The vacuum laser inten- 
sity is I - 1 x 10 19 Wcm- 2 . With n e = 1 - 10 x 10 19 
cm -3 , it gives a ratio Pl/P c = 5.3 — 53 (with Pl the 
laser power and P c = 16.2 n c /n e GW the critical power 
for self- focusing). Before focusing, the beam was split 
and 10% of the energy sent onto a BBO crystal for sec- 
ond harmonic generation. The intensity and spatial qual- 
ity of this collimated frequency-doubled beam were im- 
proved by passing through a 50-/im-pinhole-filtered beam 
reducer. It was finally incident on the jet normally to 
the main laser axis for illuminating and probing the in- 
teraction transversally (see FIG. [I]). Measurements of 
ion emission and plasma probing (interferometry and po- 
larimetry), both at 90° from the main laser propagation 
axis, were simultaneously done by means of a 45° -pierced 
mirror reflecting the probe beam onto the gas jet without 
cutting the ion beam (see FIG. HI). 

The ion detector is composed of a 100-/im spatial 
filter (solid angle acceptance is 1.2 x 10 -8 sr) with a 
Thomson parabola spectrometer coupled to 40- mm imag- 
ing microchannel plate chevron assembly (see FIG.pJa)). 
A linearized 16-bit CCD camera (1024 x 1024 pixels) 
enables single-shot recordings of scintillating traces and 
thus shot-to-shot statistical analyses of ion spectra. 

An f/1 achromatic doublet coupled with a 50/50 
beam splitter enables to record at the same time with 
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FIG. 1: Experimental setup for the interaction of a micro- 
jet and a short Ti-Sapphire pulse, illuminated by a linearly- 
polarized-frequency-doubled probe pulse, (a) Imaging Ion 
spectrometer along the normal to the main beam with Thom- 
son parabola (B = 0.3 T) and chevron microchannel plate 
(MCP). A second ion spectrometer is also set on the laser axis. 
(b) Nomarski interferometer, (c) Polarimeter using Faraday 
rotation effect. For Polarigrams R and T, the polarizers are 
rotated by respectively +10 and —10° fom the extinction 



spatial resolution ~ 1 /im, interferograms encoding the 
plasma density variations using a Nomarski interferom- 
eter [22 (see FIG. m|b)), and polarigrams revealing the 
component of the magnetic field collinear to the probe 
beam direction (see FIG.pJc)). The interferometer con- 
tains a Wollaston biprism of 10 x 10 mm 2 with a 410- 
nm interferential filter to improve the temporal coher- 
ence of the interferences. The polarimeter is composed 
of two identical assemblies, each comprising in particular 
a Glan Laser polarizer (extinction ratio better than 10 -4 ) 
and a 16-bit CCD camera. For each assembly, a scan of 
the rotating polarizer was carried out for probe beam 
transmission calibration (every five degrees over the en- 
tire rotation). For the measurements, the polarizers were 
detuned by a common angle from each minimum of trans- 
mission (typical detunings Or = +10 and 0t = —10° for 
respectively the transmitted and reflected images from 
the beam splitter, see FIG.pJc)). 

The magnetic field detection rests on the measur- 
ment of the rotation, due to Faraday effect, of the in- 
cident linear polarization of the probe beam electric 
field when it crosses the plasma [23l [24] . The total ro- 
tation angle along a path l z is given by (f) rot (x,y) = 
« — - — i n e (x.y.z)B z dz. Here B z is the magnetic field 
component collinear to the probe beam direction, n c the 
plasma critical density for the probe beam, e and m e 
the mass and charge of the electron and c the light ve- 
locity in vacuum. Taking the intensity ratio from the 
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FIG. 2: (a) and (b): MCP recordings normally to the laser 
axis for n e ~ 10 20 cm -3 , (a) short plasma gradient ~ 150 /im. 
He + trace with energy cutoff ~ 250 keV. (b) long plasma 
gradient ~ 250 /im. (c) and (d): Shadowgram at 1 ps after 
the pulse enters the jet (resolution ~ 1 /im, magnification 
8.2). (c) short plasma gradient, (d) long plasma gradient. 
The shadowgram shows a clear train of bubble-like structures 
in the laser wake. The laser polarization is normal to the 
picture (S-polarized) 

two shadowgrams R and T obtained with opposite de- 
tuning angles (see FIG.[TJc)) prevents the extracted (j) rot 
map from being corrupted by probe beam shot-to-shot 
fluctuations in intensity, and image shift from plasma 
refraction or vibrations. However, this requires a pixel- 
to-pixel correspondance of the pictures. We achieved a 
sub-pixel-accurate matching of the two pictures by tag- 
ging the field of view at three reference- areas prior to 
the beam-splitting, and by running a numerical pattern 
matching routine. 

Figure [2] shows the MCP recordings at 90° from the 
laser axis for two plasma gradient conditions. In (a), the 
laser is shot on the gas jet at an early stage of its expan- 
sion in the vacuum chamber so that the plasma gradient 
is steep ~ 150 /im (with a gaussian fit). In (b), the in- 
teraction occurs at a late time of the expansion (80 ms 
after the closing of the valve), and the gradient is smooth 
^ 250 /im. For both conditions, the peak density in the 
jet has the same value n e ~ 10 20 cm -3 . We observe 
systematically radial ion acceleration from Coulomb ex- 
plosion in case (a) and no trace in (b). In both cases, no 
ion was detected in the forward direction (see FIG. IT] for 
the geometry). When acceleration occurs, only He "are 
detected from 12 keV to cutoff energies of 250 keV. We 
attribute the absence of He 2+ , though the laser peak in- 
tensity is 1000 times higher than the helium potential of 
full ionization, to efficient recombination of this species 



in the neutral corona starting at 100/im from the laser 
axis. Indeed, the electron capture is the most probable 
channel of interaction (cross section a ~ 10 16 cm -2 ) for 
He 2+ projectile of 250 keV streaming in He gas [25]. In 
panel (d), one can see the typical shadowgram (revealing 
sharp electron density fronts) obtained in case (b). A 
train of zig-zag ordered and well separated bubble struc- 
tures with an average size of about 30 /im (^> X p = c/uo pei 
the plasma wavelength) can be clearly seen in the wake of 
the laser pulse. The bubbles have a bright spot in their 
centers that is not observable when the probe beam is 
turned off. This means that the emission is not of self- 
emission nature and results from refractive effects of the 
probe light by these structures behaving like microlenses. 
An interesting point is the quasi-circular shape of these 
structures, meaning that, during the transverse illumi- 
nation of 30 fs, their centers are quasi-static in the pic- 
ture plane. In the short-plasma scale length case, for 
which transverse ion acceleration is detected (Fig. J2ja)), 
no such isolated structure is ever observed in the gas jet 
(Fig. He)). 

To understand this inverse correlation, we investi- 
gated in time the interaction in the long gradient case. 
Figure [3] presents a sequence of five shadowgrams taken 
at different delays after the pulse enters the jet (to). Very 
early in the interaction (to + 540 fs), sizeable filaments 
develop next to the main laser channel [19] (see arrows 
in panel (b) from Fig.pl. As the pulse propagates along, 
the filamentation becomes stronger and multiple circular 
structures appear about the laser channel, near the cen- 
ter of the jet. In this time-resolved scan, we observe that 
the structures appear randomly in the laser path and no 
correspondance could be drawn between the pulse and 
bubble positions. Moreover, even if this shadowgram di- 
agnostic is not intrinsically time-resolved, contrary to the 
proton radiography technique [26 , the stability of the in- 
teraction was sufficient to assess statistically the average 
size of the structure until the moment when the shadow- 
grams become blurred and unexploitable, at about 7 ps 
after the onset of the interaction. Over that timespan, 
we observed that the structures have fairly the same size, 
irrespective of their positions to the laser front and mea- 
sured an average increase of less than 10 /im in size. 

Besides being copiously generated by the interaction, 
the nonlinear structures distribute spatially into two re- 
markable patterns, identifiable by following the bright 
central blobs, which arose randomly in our experiment. 
Figure [4] illustrates both antisymetric (zig-zag about the 
laser axis) (panel (a)) and axial (panel (b)) distributions 
of the bubbles. In Ref. [27j the authors carried out 2D 
particle-in-cell simulation with a^ = 3, r = 30 fs and 
about 1 mm of plasma at n e = 10 19 cm -3 . They ob- 
tained for initially p- and s-polarized light (i.e. an elec- 
tric field collinear and normal to the observation plane), 
respectively antisymetric and axial density map patterns 
resembling our experimental results (see figures 3 and 




FIG. 3: Time-resolved pulse penetration in the plasma gradi- 
ent, (a) reference instant to, (b) to + 540 fs, (c) to + 940 fs, 
(d) to + 1.54 ps, (e) to + 2.84 ps. Arrows indicate the growth 
of prominent laser filaments . 




FIG. 4: Spatial distribution of the bubbles for two shots in the 
same conditions at to + 1.54 ps showing (a) antisymetric (zig- 
zag) and (b) symetric axial patterns about the main plasma 
channel 

4 in [27 ). In particular, the antisymetric modulation 
is ascribed to the hosing instability (see also [9]). Al- 
though in our experiment the incident laser polarization 
is normal to the observation plane, corresponding to the 
s-polarization of simulations [9] [27] , hosing instabilities 
have also been observed experimentally in similar con- 
ditions [28]. Therefore, following [9] E3 EH], we pro- 
pose to ascribe our observations to the laser excitation 
of the hosing instability inducing the developpement of 
von Karman row of electron vortices (a) and to shadow- 
graphic traces left by post-solitons (b). 

At that point, it is clear that the interaction in 
this operating regime (n e ~ 10 20 cm -3 , long gradient 
~ 250 /im) produces numerous instabilities (filamen- 
tation, hosing, vortex/soliton) , known to seriously de- 
plete the pulse energy [6j [8] [9], and prevent efficient 
plasma channel charging for subsequent Coulomb explo- 
sion. Moreover, despite their expansions by ion core ex- 
plosion [9j[12j[T5], the bubbles do not provide a detectable 
acceleration, possibly because of a number of acceler- 
ated ions below the detector noise level. Therefore, we 
were unable to verify the post-soliton acceleration sce- 
nario predicted in [6j [15] . 



To confirm the electromagnetic nature of the coher- 
ent structures, we carried out some polametric measure- 
ments. To simplify these measurements, we have identi- 
fied a regime where few structures are generated at each 
laser shot, near the center of the jet (statistically close to 
one). For that, we used a stretched laser pulse r = 250 fs 
carrying the same energy (<2q = 0.8) and a small sonic 
nozzle of 200 jiva diameter delivering a jet of peak den- 
sity of few 10 19 cm -3 which we let expand in vacuum to 
increase the gradient before the interaction. This setup 
enables (i) a well-localized interaction, (ii) smoother in- 
terferograms and polarigrams without filaments and bub- 
ble overlapping (in) efficient energy transfer to the struc- 
tures [9]. 
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FIG. 5: (a) Shadowgram with one isolated structure in the 
laser wake at £ + 940 fs (r = 350 fs, 10 19 cm -3 ), (b) ro t map 
(in radians) of the yellow square dotted zone in (a). Detuning 
angle of the polarisers: ±10°. (c) Line-out along the black 
dotted line in (b) 

Figure [5la) shows a shadowgram with one single 
spherical structure of 57 /im diameter in a clean area. In 
that simplified case, we were able to extract a (j) rot map 
around the structure (panel (b)). From panel (b) and (c), 
we can deduce that the structure in this condition has an 
almost axisymmetric polarimetric imprint about the laser 
polarization direction. The density in the bubble shell is 
found to be ~ 7 x 10 19 cm -3 , five times that of the sur- 
rounding plasma. Though the polarimetric information 
is the convolution of B z by the density along the inte- 
grated probe beam path in the plasma, we will assume 
for the remaining that main contribution comes from the 
structure and not from the surrounding plasma, as the 
major part of the bubble seems to be lying in a weakly 
magnetized area (white region). Doing so, we thus cal- 
culate the maximum magnetic field extractable from our 
data. For n e = 7 x 10 19 cm -3 , an average path l z in 
the shell of 25 /im (from panel (c)), n c = 4 x 1.67 x 10 21 
cm -3 and a Faraday rotation of 0.04 radians, we find 



the same order of magnitude as in the numerical compu- 
tations in similar conditions [6[ [8] and is consistent with 
the fact that the fields dissipate as the structure expands. 
The extraction of the space-resolved magnetic field map 
and a refined magnetic description of the nonlinear struc- 
tures will be the subject of another publication. 

In conclusion, our reproducible observations of 
bubbles identified as nonlinear coherent structures 
(vortex/post-soliton) made it possible to demonstrate 
their detrimental effect on the efficiency of laser-ion ac- 
celaration. These structures appear abundantly in the 
laser wake when the plasma has a long gradient pro- 
file in which the laser self-focuses, filaments and develop 
hosing. When using a small nozzle with lower density, 
simple conditions with one structure per laser shot can 
be statistically obtained. We think this all-optical setup 
represents a powerful method that opens up new oppor- 
tunity of in-depth investigation of fundamental nonlin- 
ear plasma instabilities and also particle acceleration for 
users of well-spread TeraWatt-class lasers. 
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